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Abstract 

The Spectro-Radiometric Calibration Assembly (SRCA) provides on-orbit spectral 
calibration of the MODerate resolution Imaging Spectroradiometer (MODIS) reflected 
solar bands and this paper describes how it is accomplished. The SRCA has two adjacent 
exit slits: 1) Main slit and 2) Calibration slit. The output from the main slit is measured 
by a reference silicon photo-diode (SiPD) and then passes through the MODIS. The 
output from the calibration slit passes through a piece of didymium transmission glass 
and then it is measured by a calibration SiPD. The centroids of the sharp spectral peaks of 
a didymium glass are utilized as wavelength standards. After normalization using the 
reference SiPD signal to eliminate the effects of the illuminating source spectra, the 
calibration SiPD establishes the relationship between the peaks of the didymium spectra 
and the grating angle; this is accomplished through the grating equation. In the grating 
equation the monochromator parameters, (3 (half angle between the incident and 
diffractive beams) and 0 off (offset angle of the grating motor) are determined by 
matching, in a least square sense, the known centroid wavelengths of the didymium peaks 
and the calculated centroid grating angles from the calibration SiPD signals for the peaks. 
A displacement between the calibration SiPD and the reference SiPD complicates the 
signal processing. 

The response of the MODIS bands to the SRCA output is also normalized by the 
reference SiPD signal to eliminate the effect of the source spectrum. That response 
differs from what is measured by MODIS using the Spectral Measurements Assembly 
(SpMA), a laboratory double-monochromator, due to the wider slit width of the SRCA. 
The SRCA slit function, calculated using measurements by the SpMA and the SRCA at 
the same MODIS temperature, will be used to recover the MODIS spectral response 
using the SRCA data measured on-orbit. 


A series of test results, at different temperatures, indicate that the SRCA is stable and is 
capable of detecting MODIS wavelength shift on-orbit. 
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X. Introduction 


The Spectro-Radiometric Calibration Assembly (SRCA) is one of the on-board 
calibrators in MODerate resolution Imaging Spectrometer (MODIS), which is being 
designed and manufactured by Raytheon Santa Barbara Remote Sensing. The SRCA is 
capable of radiometric, spectral, and spatial calibrations. 

The reasons for on-orbit spectral calibration are: ( 1 ) Spectral band responses could 
change because of variations of the transmissive properties of the thin-film multi-layer 
filters with temperature [1]. Band shifts were observed during ground testing of 
precursor remote sensing instruments. Some on-orbit measurements have been 
interpreted as showing evident shifts in band response profiles. (2) During radiometric 
calibration of the solar reflective bands, the Sun is the main calibration source for 
MODIS. The MODIS instrument will be periodically calibrated by viewing an 
illuminated solar diffuser plate. Knowing the MODIS band shifts helps redefine the 
incoming solar radiation of the bands. (3) Band response profiles, recovered from the 
MODIS signals measured under SRCA illumination, lead to an understanding of the 
changes in the spectral response of MODIS bands. 

During spectral calibration the SRCA measures the center wavelength shifts of the 
MODIS bands. It calibrates MODIS bands 1-19 and 26 covering a wavelength range 
from visible (VIS), near Infrared (NIR), to short wavelength IR (SWIR). 


The SRCA is configured as a monochromator during spectral calibration. The SRCA 
provides end-to-end calibration because its beam passes through the complete optical 
train, but it only fills about 1/5 of the MODIS aperture. In order to ensure that the 
measurement has sufficient Signal-to-Noise Ratio (SNR), the exit slit of the SRCA 
monochromator has been widened to cover 5% to 27% of MODIS band widths 
depending upon the band. Hence, the center wavelengths measured by the SRCA are 
biased by the exit slit width effect and the partial aperture filling effect. 

However, the SRCA is expected to be a stable source. It can detect the wavelength 
change in two consecutive calibrations. During the prelaunch period MODIS is 
spectrally calibrated against a laboratory instrument — the Spectral Measurements 
Assembly (SpMA) in an ambient environment. The SpMA is a double -monochromator, 
which fills the MODIS aperture and hence provides accurate spectra for the MODIS 
bands. At the same or near the same time and environment, the SRCA is operated. The 
differences of the center wavelengths measured by the SRCA and SpMA serve as 
correction values [2]. The correction values, which are the differences between SpMA 
and the SRCA measurements, are assumed unchanged throughout MODIS life. These 
correction values are mainly due to the use of a different SRCA slit width. In the other 
measurement environments, such as in thermal vacuum testing at different temperatures 
and on-orbit, the SRCA has been and will be run repeatedly. SRCA measured center 

wavelengths compared with the previous measurements indicate band center wavelength 
shifts. 



Figure 1 is the optical layout of the SRCA, which consists of three sub-assemblies: a light 
source, a monochromator, and a folded collimator. 

The light source is a Spherical Integration Source (SIS) with six embedded lamps (four 
10W and two 1W). In the spectral mode only two lamp configurations, first three I0W 
and then one 10W, are applied to minimize the period of the spectral calibration and to 
prolong lamp life. The radiant output is monitored by a temperature-controlled silicon 
photodiode (SiPD) embedded in the SIS (not shown in Figure 1) providing a feedback 
signal to keep the lamp output constant. 

Calibration SiPD 



Figure 1 SRCA layout 

To isolate the monochromatic beam from the mix of different diffraction orders, order- 
sorting filters are inserted into the optical path. A monochromatic beam of different 
wavelengths is output when the grating rotates with the required order-sorting filters 
inserted synchronously. 

There is a secondary exit slit, which is parallel to and is located six millimeters from the 
main exit slit, with a didymium glass filter and the calibration SiPD behind it. The 
spectral beam, from the secondary slit, passes through a didymium glass and is detected 
by the calibration SiPD while the spectral beam from the main exit slit is sampled by the 
reference SiPD, with the same characteristics as the calibration SiPD, located at the 
center of the secondary mirror of the Cassegrain system. The reference and the 
calibration SiPDs are manufactured in the same lot and selected in pairs. The two sets of 
signals form the database for wavelength self-calibration. The calibration SiPD signal 
shape, after normalization by the reference SiPD signal to remove the light source 
spectral shape, results in a profile which matches the didymium transmission curve. 
However, the location difference of the two SiPDs complicates the signal processing. 



During the spectral mode heat from the lamps and motors warm the SRCA and the SRCA 
housing temperature varies with time. In order to compensate for the SiPD spectral 
sensitivity variation with temperature, the temperatures near the two SiPDs are recorded 
with the SiPD data. 


The MODIS bands spectrally calibrated by the SRCA are listed in Table 1. 


Table 1. The center wavelengths and band widths of MODIS bands 


Band 

Center wavelength (|im) 

Band width (jim) 

Band 

Center wavelength (urn) 

Band width (jim) 

1 

0.646 

0.050 

11 

0.530 

0.014 

2 

0.856 

0.045 

12 

0.547 

0.012 

3 

0.466 

0.021 

13 

0.666 

0.012 

4 

0.554 

0.021 

14 

0.677 

0.014 

5 

1.242 

0.031 

15 

0.747 

0.013 

6 

1.629 

0.036 ~ r 

16 

0.866 

0.018 

7 

2.114 

0.067 

17 

0.904 

0.040 

8 

0.412 

0.017 

18 

0.936 

0.016 

9 

0.442 

0.011 

19 

0.936 

0.055 

10 

0.487 

0.012 

26 

1.382 “1 

0.048 


2. The spectral transmi ttance of didvmium glass and peak location 

Monochromators need to be regularly re-calibrated for wavelength due to possible 
changes of instrument parameters. A line-spectrum light source is normally utilized as the 
wavelength calibrator. When the wavelength standard needs to be carried into space, a 
sharp-peaked absorption glass is a convenient choice. The didymium glass has several 
absorption and transmission peaks in the desired wavelength range. Although there is no 
guarantee that the peak location is stable, it is likely to be much more stable than the 
multi-layer thin-film interference filters used in the MODIS focal plane arrays. 


Figure 2 illustrates the spectral transmittance of the didymium glass [3]. The 
transmission profile has several sharp peaks and valleys. Selection of peaks has the 

advantage of higher SNR in 

Transmittance of didymium glass determining the centroid values 

of the peaks than selection of 
valleys. The peaks are 
numbered from D1 to D6 to 
indicate their locations. The 
sharper the peaks and the higher 
the diffraction order applied, the 
more sensitive the determination 
of the peak location by the 
calibration SiPD signal. Two 
peaks are utilized in this 
application: one at 0.496 jam 
(D2) and the other at 0.55 l|am 
Figure 2. Didymium spectral transmittance (D3). The selection of the two 

peaks considers the use of higher diffraction order and the peak locations over the ran^e 





of MODIS bands to be measured. The use of the third diffraction order for 0.496 fim 
(D32 1 ) and 0.55 lfim (D33) peaks and the second diffraction order for 0.551 pm (D23) 
locates the three peaks in the grating angle range of 6.84 to 12.06 degrees. The entire 
grating angle range for the calibration spans from 6.35 to 15.33 degrees. 

3. Spectral calibration signal and its normalization 

The basic grating equation is 

mX = A(sin 6 i - sin 6 d )cos 0 ( 1 ) 

where A is the grating spacing; m is grating order; X is the wavelength of the output 
beam; 0j=p-0 g , is the incident angle about the grating normal in the diffraction plane; and 
e d=P+0 g > is the diffraction angle about the grating normal in the diffraction plane. 0 g is 
the angle between grating normal and the bisector between incident and diffraction rays. 

(j) is the out of diffraction plane angle. 

When the grating rotates an angle of 0 M with an initial offset angle of 0 off , and assuming 
the beams to be in the diffraction plane for simplicity, the wavelength of the output beam 
from the main exit slit can be written as 

2A 

X- — sin( 0 M + 0 ')• cos p ( 2 ) 

m v 7 

if we select negative diffraction orders. 

The spectral calibration consists of determining two monochromator parameters: the half 
angle between the incident beam and the diffraction beam, (3, and the grating motor offset 
angle, 0 off , The two parameters may change under different operational environments and 
conditions. As soon as the two parameters are determined using the didymium 
wavelength standard, the relationship between grating motor angle and wavelength is 
established and then the monochromator is calibrated. 

When the monochromatic beam exits through the didymium slit, an angular distance of A 
from the main slit, Eq.(2) takes a different form: 

1 ^ sin < e M +0 off + j) • cost P + f) (3) 

In normalization, the signals of the calibration SiPD with the wavelength value of X in 
Eq. (3) should be divided by the signals of the reference SiPD with the same wavelength 
value of X in Eq. (2). In other words, the signal from the calibration SiPD at grating & 
angle of 0' M must be divided by that from the reference SiPD at grating angle of 0 M . The 
0' M corresponding to each 0 M is calculated by setting Eq.(2) equal to Eq.(3). However, 
since the exact values of (3 and 0 off are unknown before the calibration, we have to use the 
design values for the two parameters to relate 0' M to 0 M . Since the change of (3 and 0 off is 
small, the use of design values of P and 0 off introduce an uncertainty less than one step. 


D32 means a diffraction order of three and peak number two. 
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Figut-e 4. Grating angles 6' M and 0 M Figure 5. The difference between Q' M and 

in terms of wavelength q m 

Figure 4 shows the 0' M and 0 M in terms of wavelength according to Eqs. (2) and (3). The 
two curves are fairly linear and parallel to each other. The difference of 0 M - 0' M in terms 
of grating motor step is shown in Figure 5. In this application, each step is equivalent to 
0.00588 degrees on average. For spectral calibration we are interested only in the grating 
angle ranges that the didymium peaks cover. The differences (0 M - 0' M ) for the peaks 
D23 and D32 are constant (108 and 106 steps, respectively) while the peak D33 has a 
difference of mostly 106 and partly 105 steps. Considering that one step difference for 
partial didymium peak sampling does not introduce detectable error, the difference for 
peak D33 utilizes the value of 106. 


Table 2. Grating motor step number for the three didymium peaks [4] 



Didymium peaks 

^ ... 

Motor step for the 
didymium slit 

Motor step for the 
! main slit 

Number 

Wavelength 

(pm) 

Diffraction 

order 


End 



D23 

0.551 

2 

31764 

31857 

31872 

31965 

D32 

0.496 

3 

32322 

32334 

32328 

32440 

D33 

0.551 

3 

32408 

32544 

32514 

32650 


Table 2 gives the actual motor step number and the range for the three didymium peaks. 
The grating motor angle is nominally zero at the grating motor step number of 30600. 
Their differences, dstep#, are 108, 106, and 106 as mentioned. The normalized 
calibration SiPD signal, CDS Norm over the didymium peak region is 

CDS NonruL .pk(steP#,m) = ^ Nation Si PD,L,pk( Ste P#> m ) ~ calibration SiPD, dark 

DN re ference SiPD, L, pk ( step#+dstep#, m) - DN referencsSPD 

, dark 

where the subscript L is lamp configuration, ‘pk’ stands for the peak region, and m is the 
diffraction order. Dark readings are taken before the SRCA is turned on and after the 
SRCA lamps are turned off. Test data have shown that the dark reading is stable 
throughout the calibration period, an average value is assigned. 

































Figure 6 shows the signals of the calibration SiPD and the reference SiPD, both with dark 
voltage subtracted, versus grating motor step number. The locations of the three peaks 
are marked by arrows. The other reference SiPD data are applied for normalization of 
MODIS band/channel signal while the calibration SiPD data are useful only in the 
vicinity of the didymium peaks. Figure 7 shows the normalized calibration SiPD signal 
in terms of grating motor angle only in the ranges listed in Table 2. 

Before the reference SiPD signal is used to normalize the calibration SiPD signal, we 
have to handle the noise of the reference SiPD signal. The noisy reference SiPD signal is 
attributed to the high gain setting of the SiPD circuit. If the noise is not eliminated all 
calibration signals and MODIS band response profiles will be noisy after these signals 
are divided by the reference SiPD signals. A quadratic equation fits the measured 
reference SiPD to smooth out the noise without degrading the centroid value of the 
reference SiPD signal. Figure 8 shows the measured reference SiPD signal and the noise- 
smoothed reference SiPD profiles for 30W at Nominal temperature plateau (269K). 
Spurious (obviously noisy) signals have been excluded from the calibration. 
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Figure 7. Normalized calibration 
SiPD signals 


Reference SiPD data over didymium peak 
(30W at Nominal plateau) 




Figure 8. Polynomial approximation to smooth out noise of the reference SiPD signal 


4- Spectral calibratio n signal to cover all MODIS reflated solar hands 

To obtain adequate SNR levels for the 20 MODIS bands to be calibrated, two lamp 
configurations are applied: 1- 10W and 3-10W. The motor rotates the grating 






unidirectionally to avoid uncertainty due to motor backlash or motor positioning. For the 
30W lamp configuration, the grating rotation covers MODIS bands 1 - 4, 8 - 11, and 17 
plus the three didymium peaks. Then the grating rotates to its home position. For MOW 
lamp configuration, one unidirectional rotation of the grating covers wavelengths for 
MODIS bands 5-7, 12 - 16, 18, 19, and 26 plus the three didymium peaks. The two 
monochromator parameters, p and 0 off , are determined for each lamp configuration to 
account for the grating motor starting at a different offset angle when it rotates to its 
home position. The offset angle, 0 off , is critical to the calibration accuracy. 


The order sorting filter mounted on the filter wheel is positioned according to the 
diffraction order required and grating angle (Table 3). 


Table 3. The relationship of diffraction order with fill 

ter# and grating angle [2] 

Filter# 

Grating angle (degree) 

Diffraction order 

open hole 

7.8-15.0 

-1 

#1 

7.9- 15.0 

-3 

#2 

< 10.4 

-2 

#2 

> 10.4 

-3 

#3 

<8.6 

-1 1 

#3 

Note: The band pass is 0.40 - 0.52ur 

>8.6 

n for Fil tor# 1 O S'? - 0 1C\\ im fV\r 

-2 

rUl non.. c 


Filter#3. 


The recovery of 20 MODIS band responses and six didymium peaks in two 
unidirectional grating rotations requires that the grating motor start step#, ending step#, 
step size, and number of samples per collection be carefully designed. Even in this 
circumstance some bands may have to share the same grating position but need different 
diffraction orders. Hence, the order-sorting filters have to flip back and forth at the same 
grating step. 

The spectral calibration has to be executed when the earth and atmosphere are in dark to 
eliminate external-illumination effects on the high-gain reference SiPD signal. Under 
this constraint, the entire calibration has been divided into four periods; two periods per 
lamp configuration. Each period will be completed within the duration of the dark 
portion of the orbit (35.4 minutes). The operational procedure is executed according to a 
series of programmed instructions [4]. 

5. Determination of th e two monochromator parameters. [3 and 9^ 

The angular displacement of A in Eq.(3) is calculated from 

A = tan (- ) ( 5 ) 

J mono 

where dS is the distance between the main exit slit and the didymium slit (6mm), and 
fmono (260.6mm) is the focal length of the monochromator. It seems in Eq. (3) that the 
angular displacement could be another parameter so that three monochromator 
parameters are uniquely determined by the three peak positions. However, it does not 



behave that way. Notice that the parameters p and 0 off both have A/2 added to them. An 
increment of A/2 merely leads to corresponding reductions in p and 8 0ff by the same 
amount. It means A/2 can not be determined as an independent parameter. We use the 
design value of A/2. V ariation in the monochromator focal length will introduce 
uncertainty to the calibration. 

The monochromator parameters, p and 0 off , are derived by best matching the didymium 
peaks in terms of wavelength and the normalized calibration SiPD signal in terms of 
grating motor angle. The match is not the detailed profile but the centroid values of these 
peaks. Actually peak D3 (Figure 2) at 0.55 1 Jim is used with both diffraction orders 2 
and 3, these two combinations occur at two distinct grating motor angles. Hence, the two 
didymium peaks lead to three unique grating angles and also three distinct centroid 
wavelengths because of the possible difference in the number of samples which are used 
to compute the centroids for the same peak. That means we relate the three centroid 
wavelength values of the two didymium peaks to the centroid values of the normalized 
calibration SiPD signal for the three peaks through the grating equation. 


After conversion of grating motor step number to angle using the database, the 
normalized calibration signal becomes CDS nomii pk (0' M , m). The centroid value of the 
normalized SiPD signal for each peak, Q c m DS pk L (degrees), is calculated by 

Q M=Q'M2,pk,L 

]^CDS norm pk L (0 M , m)-0 M • A6 ' m 

n CDS __ d M=9'Ml,pk,L 

d m,pk,L e' M =e' M2 ,pk-L ~ ( 6 ) 

£ CDS narm pk L (Q M , m) ■ A6 m 

e M=d'Ml,pk,L 

The data points are limited by a threshold value, which equals to 0.7 CDS 
betting the threshold maintains the consistency of the calibration results. In Eq.(6), the 
0 M1 and 0' M2 are both wings of the angle values with the signal intensity greater than the 
threshold. The number of data points depends on data availability. 

The centroid wavelength of the didymium transmission peak is calculated over the 
wavelength region of X 1>pk>L and which correspond to & MhpKL and 0' M2pkL and are 
calculated by Eq.(3) with design values of P and 0 off . Sampling for the centroid 
calculation over the angle range of Eq.(6) is not necessarily equally spaced. If some 
bands have the same diffraction orders and are sampled in the peak region with different 
step size, more data points may be available and used in Eq.(6). If we ignore the situation 
and calculate the centroid wavelength of the didymium peak with a constant step size, a 
bias could be introduced and the calibration may be unstable due to the inconsistency in 
data sampling. The centroid wavelength of the didymium peak is 



(V) 


A. 


_ ‘'F 1 

'didymium, pk,L X=Aj 


X =^2,pk,L 

^ ^ didymium ( A ) ' A ’ ZiA 
pk,L 


pk,L 


where ^didyinium(A) is the spectral transmittance of the didymium glass. Because different 
lamp configurations have different data stepping strategies, the centroid wavelength of 
^'didymium, pk.L must be calculated for each individual lamp configuration. 

A least-square fit technique [5] to adjust (3 L and 0 ofa until the A^best fits A didymium>pkiL 

at a number of peaks, N peak , is in use. X c pk s L is computed by substituting Q^ D pk L into 
Eq.(3) for 0' M - 

Starting with the grating equation (3), 

m ^"didymium,pk,L ~ m ^pk,L ~ a iL ' S ^ n ^rn.pk,L + a 2L ' COS ^m,pk,L ( 8 ) 


where a, L and a 2L are the function of (3 L and 0 off L : 

ciji = 2A ■ cos(0 off L + A/2)’ cos( p L + A/2) (9) 

a 2 L “ 2A • sin( Q off L + A/2) • cos( P L + A/2) (10) 


The result of the linear fit using the three sets of Q™ s pkyL and mX didymtupkL in Eq.(8) gives 

the best predicted values of a IL and a 2L . From the Eqs.(9) and (10), the two 
monochromator parameters can be calculated. 

9 off.l= tan ~‘(—)-z) (11) 

a iL ^ 

Pl = c°s~‘( — J-) - ^ ( 12 ) 

2Acos(e ojfL + -) 2 

Figure 9 shows the match of normalized calibration signals with didymium peaks at three 
grating motor regions. The two parameters calculated at different temperature plateaus 
are listed in Table 4. The results are highly consistent. 




Figure 9. Match the centroid values of didymium peaks with that of normalized calibration 


Table 4. Monochromator parameters measured in eight calibrations 



P (half angle) 

0 otf (grating motor offset angle) 

Nominal value 

15° 

0° “1 

30W, Nominal plateau (269K) 

15.053° 

0.006° 

10W, Nominal plateau (269K) 

15.108° 

0.003° 

10W, Hot plateau (279K) 

15.059° 

0.002° 

10W, Cold plateau (253K) 

15.076° 

0.006° 

30W, ambient (1) 

15.030° 

0.000° 

10W, ambient (1) 

15.002° 

-0.001° 

30W, ambient (2) 

15.035° 

0.000° 

10W, ambient (2) 

15.045° 

0.001° 


6. Handling dark readings for suhsample data 

Three detector sizes are used in the MODIS focal plane assemblies (FPAs): 0.25km 
IFOV (Instantaneous Field of View) for bands 1 - 2; 0.5km IFOV for bands 3 - 7, and 
lkm IFOV for band numbers greater than 7. Accordingly, the sample rate varies for 
different detector sizes. The SRCA has ten frames (one frame equals to lkm EFOV) of 
data. Hence, there are four subsamples for each channel of bands 1 - 2 (for a total of 40 
samples per scan per channel), two subsamples for bands 3 - 7 (for a total of 20 samples 
per scan per channel), and one sample for band number greater than 7 (10 samples per 
scan per channel). 

The dn, digital counts after space counts is subtracted, may vary between subsamples 
because of the possible differences in subsamples such as integration time. Figure 10 
illustrates an example of the dn variation between subsamples. Since only one fully 
illuminated sample is available, this sample is picked as the signal. Notice that the dark 2 
readings vary from subsample to subsample. Sample#10 falls into the subsample#2 


The dark values in Figure 10 are those for which no light through the slit falls on the detector. 




(even sample) category, for which the dark dn is smaller than that for subsample# 1 (odd 
sample). 

The dn signal must correspond to the sample where the detector is fully illuminated by 
the slit. It appears at sample#10 in Figure 10. When the signal dn (sample#10) is picked 
up for each individual step, the dark dn is calculated by averaging over several non- 
illuminated samples with the same subsample number as the signal dn. We define it as 
dn , dn-dn^. Figure 1 1 shows dn' versus grating motor step number for that band. 3 



response using dn ' 


7,_ The center wavelengths of the MOOTS hands 

The MODIS bands respond to the SRCA monochromatic beam, which is spectrally 
shaped by the light source, grating efficiency, and the SRCA system transmittance. This 
shape must be eliminated from the measurement. The reference SiPD signal provides 
this capability. Notice that the reference SiPD signal contains not only the spectral shape 
of the lamp source but also the spectral response, R SDS (V). Hence, the normalized MODIS 
band response, R ban d,ch a nne!.NomAX must be divided by the reference SiPD signal, SDS(X), 
and multiplied by the SiPD spectral response. 

/? / 0 \ ^band. channel ( ) T) / ' ) 1 , „ ^ „ 

band. channel. Norm ' /l V SDS( X) # ^ SDS ^ (13) 

Next the subscript, channel, is defined as follows: There is a communication channel or 
data channel corresponding to each detector in each MODIS spectral band; bands 1 and 2 
each have 40 channels; bands 3-7 each have 20 channels; and bands 8-19 and 26 each 
have 10 channels. Limited by the SiPD response cut-off wavelength at about 1.1 jam, the 
SWIR bands 5 (1.242 pm), 6 ( 1.629 pm), 7 ( 2.1 15 pm), and 26 ( 1.382 pm) can not be 
normalized by the reference SiPD signal. Fortunately, this does not cause evident 
uncertainty in the calibration for these bands because: (1) the SRCA detects the center 
wavelength shift. The detection is a comparison of calibration results from prelaunch and 
on-orbit, hence, the normalization effect will be significantly reduced. (2) The 
wavelength of the SWIR is located where the lamp output profile changes slowly with 


3 11 1S antici Pate that the on-going improvement will significantly reduce the differences. However this 
algorithm will still be used. 




amp temperature. According to data analysis, a 10K change in the lamp temperature will 
cause the SRCA source profile to vary by only 0.8% and the center wavelength will 
change by only 0.03nm. 


The MODIS center wavelength is calculated by 

A V 2 dn'( band, channel, Aj , 

2 */ d Try * ^ 

X c (band, channel ) = — ^channel, Normal 

y dn'( band, channel, X) t ^ 

2-i n Try ' 

A = A i ^ band , channel. Norm ' ^ / 

where the SRCA calibration results relate the wavelength, A,, to the grating motor step 
number. 


The SRCA measured MODIS band/channel response profiles are modified by the wide 
slit width so that the spectral responses measured by the SRCA differ from those 
measured by the SpMA. The SRCA slit width is equivalent to 5% - 27% of MODIS 
bandwidths depending on the band width and the diffraction order. Figure 12 gives an 
example of band response measured by the two instruments. The SRCA measured center 
wavelength is biased, but, the SRCA measured center wavelength shift between any two 
measurements determines the MODIS band shift because the SRCA slit size is fixed and 
limited variation in system characteristics do not change the center wavelength value. 

The differences between the center wavelengths measured by the SpMA and the SRCA 
in the same environment are the SRCA wavelength correction values, dA, h H 
(£q. 10). l he MODIS center wavelengths measured over the MODIS lifetime, 

^c.band, chan nel , on-orbit > are calculated by adding the correction values to the current measured 

center wavelength, ^ c SRCA band channel on . orbit (Eq.l6), assuming that the correction values do not 
change. 
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Figure 12. Difference of measured band responses by the SRCA and SpMA 




8 , Recovery of MODIS band/ch a nnel relative spectral response profiles from SRC A 
measurements 


Recovery of MODIS band/channel relative spectral response is possible if the SRCA slit 
function does not change over time. 


The recovery methodology is to calculate the SRCA slit function by dividing the Fourier 
transform of the MODIS measured response to the SRCA by the Fourier transform of the 
MODIS measured response to the SpMA at nearly the same time and with the same 
measurement environment temperature. The light source distribution for the SRCA has 
been normalized by using the reference SiPD signal (see Eq.(13)); and the SpMA light 
source distribution has also been normalized by using measurements from a sensor of 
known response. Therefore, we can assume that the Fourier transform of the slit function 
of the double-monochromator, SpMA, is unity. The Fourier Transform of the SRCA slit 
function, S(g>), is 


C/Va ) CA ( X ) p re i aunc fi t band, channel ) 

' 'band, channel ZT/ C juf a ( 1 \ ~ 7 

t ( apMA{ ^)p re launch.band, channel ) 

where the operator, F, is the Fourier Transform. 


(17) 


This slit function (mainly depending on the slit width, grating status, stray-light level, and 
scattering) is assumed unchanged. The MODIS band/channel relative response is 
calculated using the SRCA measured profiles on-orbit by 

Band response X) onorbitbmtIchannel = F~' } ( j 8) 

( (A 'hand, channel 

Figure 13 illustrates the band response of MODIS band 8 channel 2 measured by the 
SpMA (left) and by the SRCA at Nominal temperature plateau (right) and recovered 
using the Eqs . (17) and (1 8). 



Figure 13. Spectral response recovery for band 8 channel 2 


9. Uncertainty of the spectral calibration 


Based on the measurement data, the main uncertainty sources are summarized as follows: 

(1) The uncertainty of grating spacing contains two aspects: the spacing variation with 
temperature and the spacing repeatability in manufacturing. Low expansion material 




limits the uncertainty to less than 0.2 nm for temperature change of ±20K. The spacing 
repeatability is better than 0.01 jam, which causes uncertainty of less than 0.1 nm [6]. 


(2) The grating motor angle in terms of step number shows a variation of less than a 
single step (0.00588°per step on average) [7]. The centroid calculation from more than 
29 data points significantly reduces the effect of the angular variation. 

(3) The uncertainty of monochromator parameters, (3 and 0 off , depends mainly on the 
number of data points, the noise of the two SiPDs, and the SiPD response variation with 
temperature. The number of samples over each didymium peak is between 29-35. The 
noise of the SiPD comes from possible saturation and noisy reference SiPD signal due to 
the high gain setting. The temperature correction value for the SiPD response should be 
determined based on the manufacturer’s temperature correction. 

Some MODIS centroid wavelengths can be determined for both lamp configurations 
(10W and 30W). Evaluation for this uncertainty source is from the centroid wavelength 
difference between the two calibrations at the same temperature plateau. 

(4) The focal length, f mono , varies with temperature and/or other factors, which leads to a 
variation of angular displacement, A. The uncertainty has high weight to the spectral 
uncertainty. If the focal length changes less than 1%, the uncertainty contributes the 
same order of uncertainty as the grating motor offset angle. 

(5) Reference SiPD noise is due to the high gain setting. The algorithm discards the bad 
data points and smoothes the SiPD output profiles. 

(6) SiPD spectral response changes with temperature. The temperatures are 

synchronously measured with SiPD data. Correction for this variation is part of the 
algorithm. 

(7) No normalization for MODIS signals of SWIR bands As previously 

mentioned in section 7, this uncertainty is small. 

The estimated uncertainty of the spectral calibration (precision) is given in Figure 14. It 
indicates that the uncertainty is less than the specification required for the precision in- 
flight calibration: 0.5nm for the shortest wavelength band (0.412|im) and scaled for 
wavelength greater up to ljim [8]. In Figure 14 all channels are shown for each band. 

10. Concluding remarks 

The SRCA is a unique on-board, self-calibrated, spectral source. A didymium glass 
transmission filter serves as a wavelength standard and the wavelength of the output 
monochromatic beam is scaled by determining the monochromator parameters in each 
run of the calibration. Measurements under different temperature environments indicate 
that high repeatability can be expected. The precision of the wavelength calibration 
meets the requirement of the MODIS specification. 



Uncertainty of the SRCA in spectral mode 



Figure 14. Precision of the SRCA spectral calibration 
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